Generating the curvature perturbation at the end of inflation 
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The dominant contribution to the primordial curvature perturbation may be generated at the 
end of inflation. Taking the end of inflation to be sudden, formulas are presented for the spectrum, 
spectral tilt and non-gaussianity. They are evaluated for a minimal extension of the original hybrid 
inflation model. 
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a. Introduction The primordial curvature perturba- 
tion ( exists already a few Hubble times before the 
observable Universe enters the horizon.^ Its time- 
independent value at that stage sets the initial condi- 
tion for the subsequent evolution of all perturbations in 
the scalar and adiabatic mode, which in turn seem to 
the dominant (perhaps the only) cause of inhomogene- 
ity in the Universe. The curvature perturbation is al- 
most gaussian, with an almost scale-independent spec- 
trum -Pc = (5 X 10^5)2. 

The curvature perturbation is supposed to originate 
from the vacuum fluctuations during inflation of one or 
more light scalar fields, which on each scale are promoted 
to classical perturbations around the time of horizon exit. 
Different proposals have been made regarding when and 
how this happens. In this note I explore a new possibility; 
that the significant, and possibly dominant, contribution 
to the curvature perturbation is generated at the end of 
inflation. 

b. Generating the curvature perturbation For clarity 
I assume Einstein gravity, slow-roll inflation and canon- 
ically normalized fields. None of these essential for the 
generation of the curvature perturbation. What is essen- 
tial though, to make the curvature perturbation almost 
scale- independent, is that inflation is almost exponential 
while cosmological scales leave the horizon, correspond- 
ing to a Hubble parameter H which varies slowly on the 
Hubble timescale. 

The slow-roll inflationary trajectory will lie in a sub- 
space of field space, spanned by one or more light fields 
which are the components of the inflaton. Any additional 
light fields have a negligible effect on the inflationary tra- 
jectory. 

We are going to be interested in Fourier modes of the 
curvature perturbation C with a given wavenumber k of 
cosmological interest. To evaluate them we can smooth 
the scalar fields on a scale of order 1/k (to be precise, on 
a somewhat smaller scale). This scale leaves the horizon 
at some epoch k = aH, where H = a/a is the almost- 
constant inflationary Hubble parameter. 

At horizon exit the canonically-normalized fields can 
be defined so that one of them points along the trajec- 
tory; if the inflaton has more than one component I call 



this the adiabatic component of the inflaton and denote 
it by 0, and call the orthogonal component (s) isocurva- 
ture components. If the inflaton does have more than 
one component, the direction of the trajectory changes 
slowly on the Hubble timescale. 

The vacuum fluctuation of each light field is promoted 
to a classical perturbation around the time of horizon 
exit. The perturbation 64> of the adiabatic component 
(defined just a few Hubble times after horizon exit) gives 
a contribution to C<p, which remains constant until the ap- 
proach of horizon entry. The first proposals, working in 
the context of single-component inflation, assumed that 
no further contribution to C. is generated, making the 
observed quantity 0, |2| ■ When multi-component infla- 
tion was formulated it was seen Q that the isocurvature 
components of the inflaton will generate an additional 
component to ^, but it was still assumed that ^ will have 
settled down to its observed value by the end of infla- 
tion. In other words, it was assumed that the curvature 
perturbation is still generated during inflation. 

Later, proposals were made whereby the dominant 
component of the curvature perturbation (or at least a 
significant component) is generated only at some post- 
inflationary epoch. According to most of these proposals 
the curvature perturbation is generated at, or on the run- 
up to, some reheating transition. A reheating transition 
may be defined as the decay of a more or less pressure- 
less component of the cosmic fluid (usually taken to be 
dominant) into radiation to give a practically radiation- 
dominated Universe. In the first detailed proposal Q (see 
also 0) the curvature perturbation is generated gradu- 
ally, when the initially isocurvature density perturbation 
of some 'curvaton' grows in an unperturbed radiation 
background. According to different proposals, the curva- 
ture perturbation may be generated by an inhomogeneity 
in the decay rate @ or the mass of the decaying par- 
ticle. Further proposals invoke inhomogeneous preheat- 
ing , inhomogeneous thermalization after reheating Q 
or the inhomogeneous production of topological defects 

In this paper I explore the third logical possibility, 
that a significant and perhaps dominant component of 
the curvature perturbation is generated at the transition 



^ Standard material concerning the early Universe and observation ^ The proposals of |^ 0, were first mentioned in llll . without 
is reviewed for instance in any calculation. 
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between inflation and non-inflation. 

c. General formulas During inflation, the light 
flelds (f>i correspond by definition to those directions in 
field space satisfying \riij \ <^ 1, where 



V 



(1) 



We are assuming slow-roll inflation, which is almost ex- 
ponential corresponding to H H^. The perturbations 
of the light fields can be defined on any spacetime slic- 
ing which is non-singular 0,0] in the limit H/H^ 0, 
such as a spatially flat slice. These perturbations are al- 
most massless and non-interacting. They are supposed 
to vanish at the classical level before horizon exit, cor- 
responding to the vacuum state. A few Hubble times 
after horizon exit, the vacuum fluctuation has generated 
classical perturbations (5(/)i(x), which vary slowly on the 
Hubble timescale and are practically gaussian and un- 
correlated with each other. Each of them has the same 
spectrum, Vsitu = (-fffc/27r)^, where Hk is the value of H 
at horizon exit. 

The perturbations S(j)i{'x.) determine C at horizon exit, 
and they also determine any subsequent evolution of 
Such evolution (until the approach of horizon entr y) is 
conveniently described by the SN formalism 1^ Il2l IT^ 
15j (see also ..16]). Keeping quadratic terms |l2l |. the 
time-dependent curvature perturbation, smoothed on the 
scale fc, is 



= ^7V,(fc,i)<50,(x) + i^iV,, 



iS(t)j . (2) 



Here, N(k,4>i,p) is the number of e-folds, evaluated in 
an unperturbed universe, from the epoch of horizon exit 
when the fields have specified values <j)i, to an epoch when 
the energy density has a specified value p. In the sec- 
ond line, TV, = dN/dcjy^ and TVy = d^N/d(l)id(j)j. These 
derivatives are evaluated on the unperturbed trajectory, 
on which the 4>i are determined by k making Ni and Nij 
functions of just k and the final time t. 

At some stage before nucleosynthesis, ( settles down 
to a time-independent value, which is constrained by ob- 
servation. We write down the predictions for ^ which 
follow from Eq. ||2J), as a function of t even though we are 
interested in the regime where C has settled down to its 
final value. 

Since the observed ^ is almost gaussian, one or more 
linear terms must dominate Eq. ^ giving the spectrum 



rcik,t) = J2NfikMHk/27rf. 

i 

The spectral tilt n - 1 = dlnVr/dlnk is [lllT3| 
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(3) 



(4) 



where identical indices are summed over. Here e <C 1 is 
a slow-roll parameter, defined as ^Mp{V' /V)'^ where V 
is the derivative of V along the inflationary trajectory. 
Both e and r]jm are to be evaluated at horizon exit. 

If the adiabatic component of the inflaton dominates 
Eq. |5J), the n on-g aussianity of ( is too small to ever be 
observable [13,123 • Otherwise it may be observable. The 
likely observables are the bispectrum and trispectrum, 
which alone are generated by the quadratic expansion 
(O. They are specifled respectively ^3 by quantities 
/nl and t^l. Taking the flcld perturbations to be per- 
fectly gaussian, which has been justifled for the bispec- 
trum [23,01, and ignoring the scale-dependence of V^, 
the predictions are [o.[l3l|^ 
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tnl = 2 



2 {NnNn)^ 
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(5) 
, (6) 



where A and B are of order 1 on cosmological scales. 
Present observation |22j gives roughly |/nl| < 100, and 
absent a detection the eventual bound will be [23 
|/nl| i$ 1- There is at present no bound on tnl from 
modern data, and no estimate of the bound that will 
eventually be possible. (A crude bound from COBE data 
[13 is ItnlI < 10«.) 

d. Dominance by a single field perturbation If only 
one fleld a is relevant, Eq. ^ becomes in an obvious 
notation 



leading to 



- 1 



2Va 



2e- 



Because the first term of Eq. ([7|) dominates, 
3 , _ 1 A^.. 



(7) 

(8) 
(9) 

(10) 



and ^3 "HsiL = 36/nl^/23. In this case, /nl niay equiv- 
alently be defined by writing ( = <^g — |Cg where Cg is 
gaussian."' 

If a is the adiabatic component (f> of the inflaton, 6(p is 
just a shift back and forth along the trajectory, making 



^ The general formula for tjml follows from the special cases in 
II3I but has not been written down before. For f- Nj,, the scale- 
dependenc e o f the sp ectra is taken into account in |22| . 

4 Following [il [13, llSl I am defining /nl through the curvature 
perturbation. An alternative definition of /nl works with the 
Bardeen potential. The two definitions coincide in first-order 
perturbation theory but differ at second order. 
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SN time-independent with 

M^iVl = 1 (11) 
N^^ = {2e-v)Nl, (12) 

where rj = MpV" /V is another slow-roll parameter. 
Both e and rj are to be evaluated at horizon exit in these 
expressions. 

Putting Eq. (^ij into Eqs. © and © reproduces the 
usual expressions 

n-1 = 2?7 - 6e . (14) 

Putting Eos itTTll and JTSJ into Eq. gives |/nl| < 1- 
This means "12, 20] that the non-gaussianity is too small 
to observe, as was first demonstrated in a different way 
by Maldacena ,18]. In this case the tensor fraction of 
the density perturbation is r = 16e.^ The tensor may 
eventually be observable 0| if ^ 10"''. 

If (7 is some other field, its contribution compared with 
that of <j) will be 

Then the middle term of Eq. becomes negligible giv- 
ing 

n - 1 = 2r]^a^ - 2e . (16) 

Since the tensor perturbation depends only on H the 
tensor fraction r is reduced; 

r = 16e^<16e. (17) 

It will never be observable if Vc^^/Vc^^ is bigger than 
about 10^. On the other hand, non-gaussianity given 
by Eq. ((117)) could be observable. 

e. Generating some curvature perturbation at the end 
of inflation To see how a contribution to the curvature 
perturbation may be generated at the end of inflation, 
consider first single-component inflation, where there is a 
unique inflationary trajectory V^(0), giving as a unique 
function of. proper time r up to a shift in the origin. 
All functions of (j) and its derivatives will also be unique, 
including the energy density p — V + irf^^/dr^. 

At each position inflation ends when <j) has some value 
4>c- Usually it is assumed that the end of inflation is con- 
trolled entirely by the inflaton, so that (pc is independent 



^ I adopt the current definition of r in terms of tfie spectra. An 
older definition using tlie CMB quadrupole corresponds to r = 
12.4e. 



of position. Then inflation ends on a spacetime slice of 
uniform energy density. The new possibility, considered 
here for the flrst time, is that 0o depends on some field 
cr, whose potential is so flat that it has practically no 
effect on the inflaton trajectory (consistent with the as- 
sumption of single-component inflation). Then 4)c{<j) will 
depend on position through the perturbation (5(t(x). As a 
result the change A'ba, from a spacetime slice of uniform 
density just before the end of inflation, to a spacetime 
slice of uniform density just after the end of inflation will 
have a perturbation 6N\^a = Cc(x)- (Indeed, A'ba will 
bigger, the bigger is the delay in the end of inflation.) 
This quantity C,a is the contribution to the curvature per- 
turbation generated by the end of inflation, which might 
dominate the inflaton contribution 

In the case of multi-component inflation there is a 
family of curved inflationary trajectories in, say, a two- 
dimensional fleld space. In our part of the Universe one 
of the trajectories is the unperturbed trajectory, but the 
perturbation in the fleld orthogonal to the trajectory at 
horizon exit kicks the local evolution onto nearby trajec- 
tories, causing C to vary after horizon exit. 

If C has settled down to a constant value by the end 
of inflation, corresponding to all trajectories becoming 
essentially equivalent straight lines, the generation of a 
contribution to the curvature perturbation at the end of 
inflation goes through in exactly the way that we de- 
scribed for the single-component case. In particular, 
vanishes unless some field other than the, by now essen- 
tially unique, inflaton comes into play. 

If instead the family of trajectories is still curved as 
the end of inflation is approached, the situation is quite 
different. There is now no reason for inflation to end on 
a slice of uniform density, and one expects that will be 
nonzero even if the multi-component inflaton completely 
determines the end of inflation. In the case of hybrid in- 
flation, this is because one expects different components 
of the inflaton to have different couplings to the water- 
fall fleld. In the case of non-hybrid inflation, the end of 
inflation on a given local trajectory will occur when at 
least one of the parameters e oi r] has a value of order 
1. But that value depends on the local trajectory which 
(in contrast with single-component inflation) depends on 
position, and in any case the slices of uniform e, uniform 
r] and uniform p will all be different. Again, we see that 

is expected to be nonzero, just because there is a fam- 
ily of inflationary trajectories which are still inequivalent 
as the end of inflation approaches. 

/. The sudden- end approximation To have a clean 
result, let us suppose that slow-roll inflation suddenly 
gives way to radiation domination. This case may be 
realistic in a hybrid inflation model where the waterfall 
fleld is quite heavy, and in any case the calculation should 
give a feel for the sort of effect that might be possible. We 
will denote quantities evaluated at the end of inflation by 
a subscript e. 

As slow-roll persists right up to the end, expansion 
during inflation is far more rapid than it is afterward. 
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and is practically equal to the number of extra e-folds 
of inflation. Expanding to second order in the field per- 
turbations, the perturbation in (j)e{<^) is 



cfic'Scr + ^0c"((5a-)^ . 



(18) 



During inflation, the derivatives of N{(j)) are given by 
Eqs. 1)11(1 and 1(12(1 and we denote them respectively by 
N' and N''. Then 



— Nc'(f>c'Sa 



2 

1 

+ 2 



2N,"<f>f + NU. 



Let us assume that dominates 



J' » 



(19) 
(H'-(20) 
, corresponding to 
(21) 



Then Eq. ((20(1 is dominated by the first term, giving 

2eo \ 27r 



(22) 



The spectral index is given in terms of Vaa by Eq. ((16(1 . 
The non-gaussianity is given by 1(11) (I (imposing eo ^ 1 
and |77o| <C 1) as 



NL 



(23) 



which could be observable. 

g. A concrete model Let us apply these formulas to 
a concrete model of inflation; 
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(24) 



This is the original hybrid inflation model |25j, except 
for the last term two terms. One of them specifies the 
interaction of a with the waterfall field x, and the other 
gives its potential during inflation. The interaction cor- 
responds to an inhomogcncous mass-squared for the wa- 
terfall field; 



(25) 



which is required to be negative. 

During inflation the waterfall field x is pinned at the 
origin, leaving just Kr((T) which we assume makes a light 
(|?7(T(t| ^ 1) and has negligible slope. Inflation ends when 
X is destabilized, corresponding to 



(26) 



The only effect of the extra term is to make 0c a func- 
tion of the new field a. As in the original model, we 
require 0o <C Afp, so that inflation after the observable 



Universe leaves the horizon takes place entirely in the 
regime </> ^ My. As a result e is completely negligible, 
which has two consequences. First the tensor fraction is 
unobservably small [23. Second the spectral tilt, if big 
enough to be observable, will be given by just the first 
term of the appropriate equation Q), (PJ or 1(16(1 . 
The derivatives of 0c (c) are 



A„ 



Arr ( ^ \a<y^ 



(27) 
(28) 



Putting them into Eq. gives ^o, whose spectrum is 



2ee 



27r 



The condition for C,c to dominate C^d, is^ 



A00C 



> 



(29) 



(30) 



with 77 = m^Mp/Vo- If it docs dominate, the spectral 
index is given by Eq. ^ and has nothing to do with 77. 
Therefore 77 can have any value ^ 1, and for a crude 
order of magnitude estimate it should even be possible 
to consider 77 ~ 1. As a result the right hand side of 
Eq. ((30(1 can be far below 1, and so can the left hand 
side. 

Let us assume that dominates. Then the spectral 
tilt, if observable, is ri — 1 = ^rjaa- The non-gaussianity 
parameter, if observable, is 



5^^^ - 2 {-Ka^ 



(31) 



The predictions of the model involve the unperturbed 
field value a and its effective mass parameter rj^^. The 
unperturbed value should be taken as the mean of cr(x) 
within the observable Universe, as opposed to its value 
in some exponentially larger region |l3|. This situation 
has been discussed previously for the curvaton 0,^3 (see 
also 26] for the axion) . As cr is light it will not be trapped 
at a minimum of the potential even if inflation has lasted 
for a long time before the observable Universe enters the 
horizon. However in that case there will be a calculable 
probability distribution for a, and one is then free to as- 
sume that the observable Universe is typical which gives 
the rough value of a. Because inflation does not 'know' 
about a (until the very end) it may be regarded as quite 
likely that the lightness condition \rjcrcr \ <C 1 will be very 



® Remember that A^e, the value of N at the end of inflation, is 
called simply A'^ when one discusses unperturbed inflation with 
a value A'^ ~ 60. I am assuming V ~ Vb as in the original model. 
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well satisfied, making the spectral tilt too small to ob- 
serve. (This is in contrast with many inflation models 
where r] increases as inflation proceeds giving a spectral 
tilt of order 1/N H^]-) If the spectral tilt is big enough 
to observe there is some prejudice towards n > 1 because 
the peak of the probability distribution is at the mini- 
mum of V{a), but n < 1 is also quite reasonable at least 
if (T is a pseudo-Goldstone-Nambu-boson corresponding 
to a periodic potential. 

h. Conclusion The existence of many scalar fields is 
mandatory if the field theory is supersymmetric, and is 
anyhow suggested by string theory even if the field the- 
ory is non-supersymmetric. Flat directions in field space, 
corresponding to a small quadratic term in the potential 
and a negligible quartic term, can occur in supersym- 
metric and conformal 28] field theories. As a result, it 
is reasonable to suppose that the vacuum fluctuation of 
several scalar fields could be promoted at horizon exit to 
become a classical perturbation, and to expect that one 
will find a variety of possible mechanisms for generating 
the curvature perturbation. 

It has been known for a long time that the curvature 
perturbation changes with time if and only if the func- 
tion P{p) is inhomogeneous. It follows that a (more or 



less sharp) phase transition, corresponding to a change in 
this function, will generate some curvature perturbation 
if and only if it occurs on a spacetime slice which is not 
one of uniform p. That is one possible way of generating 
some of the curvature perturbation, or practically all of 
it. The well-known example @ is a reheating transition 
(P ~ to P = p/S), and I have explored another exam- 
ple which is the transition at the end of inflation (P ~ —p 
to P < p/3). Judging by the hybrid inflation example, 
the new mechanism is very easy to implement, and can 
take place for a wide range of the parameters. 

After the first version of this note was released, I was 
made aware that a different model generating a contri- 
bution to the curvature perturbation at the end of hy- 
brid inflation was proposed in The model invokes 
an inhomogeneous coupling instead of an inhomogeneous 
mass. Also, the calculation is done by matching metric 
perturbations whereas I am using the 6N formalism. 
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